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A  NEW  MODEL  OF  THE  GEOMAGNETIC  VARIATION  IN  THE  UPPER  ATMOSPHERE 


Jack  W.  Slowey 

Harvard-Smithsonian  Center  for  Aatropbysics 
Cambridge,  Massachusetts 


1.0  INTRODUCTION 


Large  variations  in  the  density  of  the  thermosphere  and  exosphere  in 
association  with  geomagnetic  disturbances  were  first  deteoted  by  Jacchia 
(1939,  1931)  in  the  drag  on  artificial  satellites.  Densities  obtained  from 
the  analysis  of  drag  were  then  used  extensively  to  study  and  model  this 
geomagnetlo  variation  in  the  upper  atmosphere.  Model  development  based  on 
drag  analysis  culminated  with  the  model  derived  by  Jacobia  fit  fil  (1967) 
that  was  subsequently  incorporated,  with  only  slight  modification,  in  the 
1972  CIRA  reference  atmosphere  and  other  general  models  of  the  upper  atmos¬ 
phere  derived  from  drag  analysis. 

The  models  of  the  geomagnetlo  variation  based  on  drag  analysis  were 
extremely  simple  in  form  and  usually  represented  the  variation  as  being 
uniform  over  the  globe.  As  it  turns  out,  the  geomagnetlo  variation  is 
extremely  complex  in  form  but,  while  this  was  oertainly  suspected,  the  drag 
data  revealed  few  details  beoause  of  their  poor  resolution.  It  was  not 
until  data  from  1q  situ  accelerometer  and  mass  speotrometer  experiments 
became  available  that  the  geomagnetic  variation  could  be  modeled  in  any 
detail.  The  process  is  still  far  from  being  complete,  however. 

There  are,  at  present,  two  models  based  on  high-resolution  data  in 
use.  One  of  these  is  that  of  Jaoohia  fit  fil  (1976,  1977)  that  is  incorpo¬ 
rated  in  Jacchia 's  most  recent  general  model  of  the  thermo  where  and  exo¬ 
sphere  (1977) .  The  other  is  that  of  Hedin  fit  fil  that  is  incorporated  in 
the  MSIS  model  (1977a,  1977b,  1979).  Here,  we  will  describe  initial  ef¬ 
forts  directed  toward  extending  and  improving  the  model  of  Jaoohia  fit  fil. 
In  this,  we  have  four  main  aims.  These  are:  * 

f  1)  To  incorporate  the  variation  with  looal  magnet io  time.  As  it  is 

now,  the  model  lnoludes  only  the  mean  variation  with  latitude. 
The  looal  time  variation  is  extremely  complex,  but  the  main  fea¬ 
tures  are  quite  significant  and  oan  be  inoluded  without  too  muoh 
difficulty. 
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2)  To  allow  the  form  of  the  variation  to  vary  with  the  level  of 
disturbance.  Slowey  (1981)  has  shown  the  latitudinal  forte  of  the 
variation  to  broaden  in  high  latitudes  with  increasing  levels  of 
disturbance.  This  effect  appears  to  be  related  to  a  shift  in  the 
location  of  soae  of  the  heat  input.  Mecg  (1982)  has  reported  that 
the  polar  cusp  region  shifts  significantly  toward  lower  latitudes 
for  larger  disturbances. 

3)  To  allow  for  the  prior  heat  input  in  establishing  an  index  for  the 
level  of  disturbance.  The  Kp  geomagnetic  index  admittedly  leaves 
much  to  be  desired  as  an  indicator  of  disturbance  in  the  atmos¬ 
phere.  It  is  increasingly  dear,  however,  that  there  is  much  to 
be  gained  by  taking  acoount  of  the  persistence  of  the  effects  of 
disturbance.  Theoretical  calculations  such  as  those  by  Fuller- 
Rowell  and  Rees  (1981)  and  the  semiempirioal  models  of  Hedin  st.  fll 
(1981)  indicate  that  the  disturbance  history  over  the  preceding  S 
days  is  probably  significant. 

4)  To  correct  what  is  now  seen  as  a  rather  serious  flaw  in  the 
original  model  in  not  permitting  an  inorease  in  exospheric  temper¬ 
ature  at  the  equator.  The  effect  at  the  equator  was  modeled  en¬ 
tirely  In  terms  of  an  "equatorial  wave”  of  density  that  was  as¬ 
sumed  to  originate  in  high  latitudes  and  progress  towards  the 
equator.  While  it  is  clear  from  the  mass  speotrometer  data  that 
suoh  a  wave  is  a  significant  component  of  the  low  latitude 
response,  it  is  now  equally  clear  that  a  thermal  component  also 
exists  in  low  latitudes.  Separating  the  tiro  components  in  data 
from  a  single  height  or  from  a  limited  range  of  height,  whioh  was 
all  that  was  available  to  us  until  recently,  is  extremely  diffi¬ 
cult,  however,  and  a  final  determination  of  the  relative  impor¬ 
tance  and  extent  of  the  equatorial  wave  will  depend  on  future 
analysis  of  data  from  lower  heights. 


2.0  NATURE  OF  THE  VARIATION 


In  figure  1  we  show  a  plot  of  molecular  nitrogen  density  measured 
during  a  portion  of  one  orbit  of  the  Atmosphere  Explorer-C  satellite  by  the 
OSS  aass-apeetrometer.  The  data  were  oolleoted  during  a  period  of  rela¬ 
tively  high  geomagnetic  aotivlty  (1.5  hours  earlier,  Kp  was  6+)  and  when 
the  satellite  was  in  a  nearly  circular  orbit.  Sinoe  the  height  varied  so 
little,  in  was  safe  to  reduoe  them  to  a  single  height  for  purposes  of 
analysis  and  that  is  the  way  they  are  shown  here.  The  abscissa  is  the  time 
in  minutes,  starting  from  an  arbitrary  zero  point.  The  scale  of  the  or¬ 
dinate  of  the  Ns  number  density  is  not  shown,  but  it  represents  the  ohange 
in  the  logarithm  of  the  number  density.  The  maximum  density  here  corre¬ 
sponds  to  a  ohange  of  0.5  in  log*,  n(Ns),  an  Inorease  by  a  faotor  of  more 
than  3  in  number  density,  relative  to  the  density  at  the  origin  of  the 
plot.  The  smooth  ourve  at  the  top  of  the  figure  is  the  (adiabatic  invari¬ 
ant)  geomagnetlo  latitude  plotted  on  the  scale  at  the  left  of  the  figure. 
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Figure  1.  Logarithm  of  Ns  density  at  •  height  of  250  la  as 
■assured  fay  the  OSS  mss  speotroaeter  on  AE~C  during  a  portion  of 
orbit  4919.  The  ebsolssa  Is  the  tlae  in  minutes  from  an  arbi¬ 
trary  zero  point  and  the  ordinate  la  not  given  (see  text).  The 
saooth  curve  is  the  invariant  geomagnetic  latitude  with  ordinate 
soale  on  the  left. 


The  strong  dependenoe  of  the  geoaagnetio  variation  on  geoaagnatio  lat¬ 
itude  is  well  established  and  we  can  see  It  in  orbit  after  orbit  of  the  AE- 
C  and  other  satellite-borne  aass  speotroaeter  data*  not  only  In  N,  but  in 
all  constituents.  Not  all  constituents  behave  the  way  M,  does*  however* 
and  this  is  another  iaportant  aspeot  of  the  geoaagnetio  variation.  The 
|  behavior  of  the  lighter  constituents*  suoh  as  He*  is  just  the  opposite  of 

that  of  N,  j  their  nuaber  densities  deorease  where  those  of  Nt  lnoreaae. 
And  the  density  of  atoalo  oxygen*  which  Is  of  primary  laportanoe  as  far  as 
total  density  is  oonoerned  throughout  aost  of  the  thermosphere*  aay  either 
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Increase,  decrease  or  remain  unchanged.  Zn  addition  to  a  purely  thernal 
effect,  there  is  another  factor  involved  in  whioh  the  aoleoular  mass  is 
very  important.  We  shall  represent  this  faotor  in  terms  of  a  variation  in 
the  height  of  the  homopause. 

The  variation  is  also  dearly  not  just  dependent  on  the  geomagnetic 
latitude.  There  are  two  distinct  maxima  in  the  data  of  figure  1  that 
represent  separate  regions  of  enhancement  in  the  vicinity  of  the  auroral 
oval.  As  the  figure  demonstrates,  there  is  also  an  appreciable  over-all 
asymmetry  with  respeot  to  geomagnetic  latitude. 


3.0  FORM  OF  THE  MODEL 


The  model  of  Jacchia  si  hi  represents  the  density  and  composition 
changes  that  occur  in  association  with  a  geomagnetic  disturbance  by  a  local 
increase  in  exospheric  temperature  and  a  proportional  increase  in  the 
height  of  the  homopause .  The  increase  in  the  height  of  the  homopause  is  a 
convenient  device  by  which  to  represent  an  effect  (the  inverse  variation  of 
the  lighter  constituents)  that  may  actually  be  due  more  to  wind-induced 
vertical  diffusion.  Superimposed  on  these  two  effects  is  the  "equatorial 
wave"  referred  to  in  the  introduction,  in  which  the  number  densities  of  all 
constituents  inorease  in  the  same  proportion  and  that  is  oentered  on  the 
equator.  We  will  adopt  the  same  representation  here  and  model  the  change 
in  the  logarithm  of  the  number  density  of  the  speoies  i  as  the  sum  of  three 
components: 


Aglngn^  -  Axlogni+Aglogni+Aelogni  (1) 

where  Af  log  n^  is  the  thermal  component,  Ag  log  a*  is  the  component  due  to 
the  change  in  the  height  of  the  homopause  and  Ae  log  n^  is  the  component 
due  to  the  equatorial  wave.  The  thermal  component  Aj  log  n^  is  to  be 
evaluated  from  an  atmospheric  model  assuming  an  increase  in  exospheric  tem¬ 
perature  AqT.  given  by 


AGT.-AF<0,X)  (2) 

where  A  is  the  amplitude  given  by 

A  ■  37.5Kp[l  +  0.027  exp(0.4Kp)l  ,  (3) 

as  in  the  earlier  model,  and  F(0,X)  is  given  by 
F(0»  X)  -  aca  +a,»  sinn0  + 

cos*  0(aaa  sin  X  +  axl  oosX)  + 
sin20(ala  sin  X  +  aaa  cosX)  + 

sin4  0  sin  20 (aaa  sin  X  +  aa a  oos X) ,  (4) 


4 


p^1 


t 


■» 


V 


* 


where  0  and  X  are  the  geomagnetic  latitude  and  looal  magnetic  time,  respec¬ 
tively.  To  take  the  effects  of  persistence  into  account,  we  assume  Kp'  in 
equation  3  to  be  the  weighted  mean  of  the  lagged  3-hourly  Kp  geonagnetio 
index  taken  over  the  41  values  in  the  preceding  S-day  interval  as  follows: 


Kp(t) 


where  t  is  the  time  in  days. 


t5Kp(ti)e“c(t“ti“T) 

,idh l5 _ 

t-c  y 

^  “0 (t-ti~T) 
t^“t— c-5 

c  -  1.0  d-1  and  x  is  the  time  lag  given  by 


(5) 


t  -  0.05  +0.1  eos*  0  (day) 


(6> 


In  equation  4,  the  first  two  terms  represent  the  mean  latitudinal  variation 
of  the  Increase  in  exospheric  temperature.  The  remaining  three  terms 
represent  the  variation  of  the  temperature  increase  with  looal  magnetic 
time  (LMT)  in  low,  middle  and  high  latitudes,  respectively.  To  Introduce  a 
variation  in  the  mean  shape  of  the  latitudinal  variation  with  the  level  of 
disturbance,  we  take  n  in  equation  4  as 


n  ■  5.0-Kp/3.0  (7) 

A  preliminary  determination  of  the  remaining  parameters  in  equation  4  was 
made  by  a  least  squares  fit  to  values  of  the  exospherio  temperature  in¬ 
crease  inferred  from  measurements  made  by  the  ESR04  mass  spectrometer. 
Exospheric  temperatures  were  obtained  from  N,  densities  by  Inverse  interpo¬ 
lation  in  an  atmospfaerlo  model.  The  temperature  inorease  was  taken  to  be 
the  difference  between  the  exospherio  temperature  corresponding  to  the  ob¬ 
served  Na  density  and  that  corresponding  to  the  density  computed  for  quiet 
conditions  from  the  ESR04  "quiet  time"  model  of  von  Zahn  fit  (1977).  The 
resulting  coefficients  are 
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The  isotherms  of  the  relative  temperature  inorease  as  given  by  equa¬ 
tion  4  in  the  case  where  Ki  -  3  (n  **  4)  are  plotted  over  the  globe  in 
geomagnetic  coordinates  in  figure  2.  As  oan  be  seen,  there  is  a  considera¬ 
ble  asymmetry  with  respect  to  the  pole.  The  maximum  is  between  6b  and  9h 
at  a  latitude  of  about  80  degrees.  This  probably  refleets  both  Joule 
heating  by  the  westward  auroral  electrojet  and  partiole  precipitation  in 
the  ousp  region.  The  extent  of  the  asymmetry  oan  perhaps  be  better  seen  in 
figure  3,  where  the  profiles  for  6b  and  18b  LMT  -  near  the  extremes  -  are 
plotted  together  with  the  mean  latitudinal  profile.  Of  oourse,  the  differ¬ 
ence  in  high  latitudes  as  seen  by  an  orbiting  objeot  would  tend  to  be 
smoothed  out  over  intervals  on  the  order  of  a  day  beoause  of  the  earth's 
rotation.  This  would  not  be  the  oase  in  low  latitudes,  however,  and  it  is 
interesting  to  note  that,  even  with  the  poor  resolution  of  satellite  drag. 
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Roemer  (1971)  was  able  to  deteot  this  asymmetry  in  densities  deduced  from 
drag.  -  - - 
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Figure  2.  Global  isotherms  of  the  relative  temperature  lnorease 
due  to  geamagnetlo  disturbance  as  given  by  equation  4  for  Kp  ■>  3. 
The  coordinates  are  geomagnetic  latitude  and  magnetic  local  time. 


We  should  point  out  that  a  geomagnetic  disturbance  will  alter  the 
atmospheric  temperature  profiles,  so  that  the  procedure  we  have  outlined  of 
entering  exospheric  temperatures  in  a  static  model  to  determine  the  effects 
on  the  number  densities  is  not  striotly  valid.  There  will  be  a  height 
dependence  that  we  have  not  accounted  for  that  beoomes  an  important  consid¬ 
eration  at  lower  heights.  Zn  this  connection,  we  have  developed  disturbed 
temperature  profiles  that  give  good  results  in  representing  observed  den¬ 
sity  variations  at  heights  as  low  as  150  km.  We  have  also  developed  ex¬ 
pressions  to  represent  the  effects  of  these  disturbed  profiles  on  the  num¬ 
ber  densities  analytically.  These  expressions  are  rather  involved,  so  we 
will  not  repeat  them  here,  but  they  represent  an  integral  part  of  the 
thermal  component  of  the  geomagnetic  variation  as  we  would  model  it. 

Concerning  the  other  two  components  of  the  geomagnetic  variation,  we 
would  include  them  in  the  form  given  recently  by  Slowey  (1983).  Specifi¬ 
cally,  we  would  compute  the  component  due  to  the  change  in  the  height  of 
the  homopause  from 


Ag  log  n*  «  a^AZ{], 


(9) 
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Figure  3.  Latitude  profiles  of  the  relative  temperature  increase 
for  Kj,  -  3  as  given  by  equation  4  for  6h  and  18h  LMT.  The 
profile  of  the  mean  latitudinal  variation  is  also  shown. 
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where  Azg  (meters)  Is  to  be  computed  from 

Azjj  -  22.0  AqT.  , 

where  Aqt-  is  given  by  equation  2,  and  the  oj/s  are: 

«(Ar)  *  +3.07  x  10“*  (mks) 
u(0,) -+1.03x10“*  (mks) 
o(N#)  -  0.0 

a(0)  --5.75x10“*  (mks) 
a(He)  -  - 6.30x  10-*  (mks). 

The  component  due  to  the  equatorial  wave  would  be  computed  from 


(10) 


(11) 


log n^  -  ^ log  p  ■  5.2x10“  *  loos*  t  (12) 

where  p  is  the  total  density  and  i  is  given  by  equation  3.  As  was  men¬ 
tioned  in  the  introduction,  some  work  remains  to  be  done  to  better  separate 
this  component  from  the  thermal  component  in  low  latitudes.  Suitable  data 
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are  available  from  both  mass  spectrometer  and  accelerometer  experiments. 
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